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Предложены новые способы расчета фоновых концентраций веществ в поверхностных водах и нормативов их предельно допу- 
стимых сбросов в реки, озера, водохранилища и болота со сточными водами. Создана математическая модель фонового хими- 
ческого состава поверхностных вод. Разработана методика нормирования сбросов сточных вод в водные объекты. 


Введение 

В соответствии с Водным кодексом Российской 
Федерации, целями государственного управления 
водными ресурсами являются: 1) обеспечение прав 
граждан на чистую воду и благоприятную водную 
среду; 2) поддержание оптимальных условий водо- 
пользования; 3) сохранение качества природных 
вод в состоянии, отвечающем санитарным и эколо- 
гическим требованиям; 4) защита водных объектов 
от загрязнения, засорения и истощения; 5) предот- 
вращение или ликвидация вредного воздействия 
вод; 6) сохранение биологического разнообразия 
водных экосистем. Основным способом, обеспечи- 
вающим достижение указанных выше целей, явля- 
ется нормирование допустимых воздействий хо- 
зяйственной деятельности на поверхностные вод- 
ные объекты. Практическая реализация этого спо- 
соба возможна лишь на основе разработки норма- 
тивов предельно допустимых сбросов (ПДС) за- 
грязняющих веществ в поверхностные водные 
объекты, призванных смягчить последствия одного 
из видов водопользования - отведения сточных 
вод (водоотведения). 

Согласно [1], «нормативы ПДС устанавливают- 
ся для каждого выпуска сточных вод проектируе- 
мых (реконструируемых) и действующих пред- 
приятий-водопользователей, исходя из условий 
недопустимости превышения предельно допусти- 
мых концентраций вредных веществ (ПДК) в кон- 
трольном створе или на участке водного объекта с 
учетом его целевого использования, а при превы- 
шении ПДК в контрольном створе — исходя из 
условия сохранения (не ухудшения) состава и 
свойств воды в водных объектах, сформировав- 
шихся под влиянием природных факторов. При 
этом учитываются ассимилирующая способность 
водного объекта и оптимальное распределение 
массы сбрасываемых веществ между водопользова- 
телями, сбрасывающими сточные воды». В соот- 
ветствии с тем же документом, величина ПДС от- 
дельного выпуска сточных вод есть произведение 
их расхода д т на допустимую концентрацию за- 
грязняющего вещества С дк (ПДС =< 7 ш С да ). Обычно 
предполагается, что величина д т задана оптималь- 
но, исходя из действующих нормативов водоотве- 
дения. Поэтому основной задачей определения 
нормативов ПДС является объективная оценка до- 


пустимой концентрации веществ в сточных водах. 
Однако существующие способы расчета С ДК не бе- 
зупречны как с теоретической, так и с практиче- 
ской точек зрения, что определяет актуальность их 
дальнейшего совершенствования с учетом опыта 
нормирования антропогенных воздействий на тер- 
ритории Российской Федерации в течение послед- 
них двух десятилетий. 


Существующие способы расчета ПДС 


Используемые в настоящее время способы рас- 
чета нормативов ПДС базируются на упрощении 
системы дифференциальных уравнений Навье- 
Стокса, неразрывности потока жидкости, тепло- и 
массопереноса в водной среде и донных отложе- 
ниях, причем уравнение массопереноса записыва- 
ется отдельно для растворенной, взвешенной и 
коллоидной форм миграции элементов [2]. Полу- 
чить решение указанной системы уравнений - за- 
дача исключительно сложная, в связи с чем в ряде 
случаев рассматриваются система уравнений 
Сен-Венана, одномерные уравнения диффузии ве- 
щества в водном потоке и донных отложениях. Но 
использование и этой модели сопряжено с много- 
численными затруднениями, включая недостаточ- 
ное количество необходимой исходной информа- 
ции [3). Поэтому с целью получения приемлемых 
для практического использования расчетных зави- 
симостей обычно рассматривается уравнение мате- 
риального баланса, преобразуемое к виду: 


С Р 


ф Оф / ”Е СстС/ст 

(2фУ т ср . т 


( 1 ) 


где у - коэффициент смешения поверхностных и 
сточных вод; С ф и 0 Ф - концентрация вещества в во- 
де и расход воды в водном объекте до поступления 
сточных вод с расходом д ст и содержанием вещества 
С от ; С р - результирующая концентрация вещества в 
водном объекте в контрольном створе, удаленном от 
выпуска стоков не далее, чем на 500 м [4, 5]. 

Полагая, что деструкция вещества в водном по- 
токе описывается линейным дифференциальным 
уравнением первого порядка, производится совме- 
щение решения этого уравнения и уравнения (1), 
результатом чего является выражение (2): 

Ср = Сф @* у + Сст Ч ст еХ р {-кі), (2) 

С'.фУ “I - Цст 
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где 1 - координата времени; к - удельная скорость 
изменения концентрации вещества (для так назы- 
ваемых консервативных веществ обычно принима- 
ется к= 0). Ур. (2) после замены С)=ПДК и 
п=( ОфУ~^Чт)/Чт нетрудно преобразовать к более при- 
вычной форме (3), используемой для расчета допу- 
стимых концентраций вещества в сточных водах С ДК : 

Сдк = и(ПДК ехр(/сД -Сф) + Сф. (3) 


Параметр п - кратность разбавления - в ур. (3) 
может быть получен методом Фролова-Родзиллера 
или по выражению: 


п = 


С с т Сф 
Смаке Сф 


(4) 


где С шкс - максимальная концентрация вещества в 
контрольном створе, определяемая по методу 
А. В. Караушева [4] или иным способом решения 
уравнения диффузии в водном объекте. Описанная 
выше модель, на первый взгляд, отличается просто- 
той и характеризует основные процессы антропо- 
генной трансформации химического состава по- 
верхностных вод за счет их смешения со стоками и 
деструкции загрязняющего вещества. Однако эти 
положения не совсем верны из-за целого ряда при- 
чин. Перечислим лишь некоторые из них. Во-пер- 
вых, простота модели исчезает при попытке учесть с 
помощью многочисленных коэффициентов разно- 
образные природно-техногенные условия, наблю- 
даемые на реальных водных объектах. Во-вторых, 
изменение концентраций веществ далеко не всегда 
может быть описано указанным выше способом. 
В-третьих, в настоящее время отсутствует теорети- 
ческое обоснование способа расчета параметра С ф в 
ур. (3), имеющего смысл фоновой концентрации. 

Согласно [1], фоновой называется концентра- 
ция, «рассчитываемая применительно к данному 
источнику примесей в фоновом створе водного 
объекта при расчетных гидрологических условиях, 
учитывающая влияние всех источников примесей 
за исключением данного источника». На практике 
за фоновую концентрацию вещества С ф принима- 
ется статистически обоснованная верхняя довери- 
тельная граница среднего арифметического содер- 
жания этого вещества, рассчитанная с доверитель- 
ной вероятностью 95 % по результатам гидрохими- 
ческих наблюдений в створе, расположенном вы- 
ше по течению от выпуска сточных вод. Если меж- 
ду расходом воды и концентрацией рассматривае- 
мого вещества существует статистически значимая 
связь, то расчет проводится по формуле (5), а если 
нет - то по формуле (6): 

С*=С(0,5 % ) + ^|р (5) 

Сф=С м +Ц^, ( 6 ) 

ѵм 

где С(0 95% ) - концентрация вещества, определенная 
по регрессионной зависимости при минимальном 
среднемесячном расходе воды в год с водным стоком 
обеспеченностью 95 %; І ка — критическое значение 


коэффициента Стьюдента при уровне значимости 
а = 5 %; М— длина ряда гидрохимических наблюде- 
ний; а - среднее квадратическое отклонение кон- 
центрации вещества; У - средняя квадратическая 
погрешность расчета концентрации вещества от рас- 
хода воды по регрессионной зависимости; С м - сред- 
нее арифметическое значение концентрации веще- 
ства в водном объекте в период года, наиболее небла- 
гоприятный с точки зрения самоочищения [6]. 

Процедура вычисления фоновых концентраций, 
как и расчет параметра п, существенно осложняется 
использованием ряда поправок и допущений, в ре- 
зультате чего на практике их определение нередко сво- 
дится к использованию формулы (6) для зимнего пе- 
риода. Но самое главное, в настоящее время нет мето- 
дики определения антропогенной и природной соста- 
вляющих С ф , особенно для случаев, когда фоновая 
концентрация больше предельно допустимой. По этой 
причине громоздкие вычисления оказываются беспо- 
лезными, так как в нормативных документах указыва- 
ется на целесообразность в подобных случаях прово- 
дить расчет нормативов ПДС на основе предельно до- 
пустимых концентраций ( С да =ПДК). Кроме того, сле- 
дует отметить, что в настоящее время разработан и 
действует на всей территории РФ нормативный доку- 
мент только для расчета фоновых концентраций в во- 
дотоках, методики расчета фоновых концентраций в 
озерах и водохранилищах отсутствуют, а расчет фоно- 
вых концентраций в болотных водах регламентирован 
только на территориях отдельных областей. 

С учетом актуальности перечисленных вопро- 
сов, неоднократно обсуждавшихся в работах 
Е.В. Веницианова, А.К. Кочаряна, Д.Я. Ратковича 
и других авторов, на кафедре гидрогеологии, инже- 
нерной геологии и гидрогеоэкологии (ГИГЭ) Том- 
ского политехнического университета (ТПУ), в 
Областном государственном унитарном предприя- 
тии «Территориальный центр Томскгеомонито- 
ринг» (ТЦ Томскгеомониторинг), Областном госу- 
дарственном учреждении «Облкомприрода» 
(г. Томск) и Томском отделе водных ресурсов Верх- 
не-Обского бассейнового водного управления на- 
чаты исследования с целью разработки математи- 
ческих моделей природно-антропогенной тран- 
сформации состояния водных объектов, пригод- 
ных для решения задач расчета нормативов ПДС и 
оценки воздействия на окружающую среду в запад- 
но-сибирском регионе. В рамках этих исследова- 
ний при непосредственном участии автора были 
разработаны «Временные методические указания 
по проведению расчетов фоновых концентраций 
веществ в болотных водах и предельно допустимых 
сбросов (ПДС) вредных веществ в болота со сточ- 
ными водами», введенные в действие на террито- 
рии Томской области с 23.07.2003 г. В данной рабо- 
те приводятся результаты следующего этапа иссле- 
дований, связанного с попыткой унификации ме- 
тодик расчета нормативов ПДС для всех типов по- 
верхностных водных объектов, разработки и обос- 
нования нового способа определения фоновых 
концентраций веществ в поверхностных водах. 
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Оценка допустимых антропогенных воздействий 

Расчет допустимой концентрации вещества С дк в 
сточных водах проводится для случая наиболее не- 
благоприятных гидрологических условий, когда на- 
блюдается минимальное самоочищение поверх- 
ностных вод от загрязняющих веществ. Но при этом 
на практике не учитываются погрешности опреде- 
ления исходных параметров ур. (3) и однозначно не 
определены действия, когда С ф больше ПДК. Не ме- 
няя формулировок нормативных документов, выход 
из указанной ситуации может быть найден при ис- 
пользовании для определения С дк традиционных 
критериев однородности двух выборок объемом 
М— в условно фоновом (Сф) и нарушенном (ф со- 
стояниях. Возможность подобного подхода связана 
с тем, что смысл расчета нормативов ПДС как раз и 
состоит в нахождении такого уровня антропогенно- 
го воздействия на водный объект, при котором его 
состояние существенно не меняется. Иными слова- 
ми, вероятность Р отвергнуть нулевую гипотезу # 0 
об однородности двух выборок не должна превы- 
шать принятый уровень значимости а: 

Р(Е, е В | Но) <а, 

где Е, - статистика для проверки нулевой гипотезы; 
В- критическая область. Если предположить, что зна- 
чения С ф и С р распределены по закону Гаусса, а диспер- 
сии (Уф и <у] известны и равны (а 2 =а ф =аІ при постоян- 
ном водоотведении), то в качестве статистики возмож- 
но использование критерия Стьюдента I в виде: 



или с учетом (1): 



Тогда ограничение для величины С дк при задан- 
ном критическом значении 2 ка может быть сфор- 
мулировано в виде: 

Сдк <Сф + п2к г аСТ ( 7 ) 

Таким образом, допустимая концентрация ве- 
щества в сточных водах с поправкой на кратность 
разбавления п не должна выходить за пределы 
верхней доверительной границы изменения фоно- 
вого содержания в поверхностных водах. С учетом 
этого задача сводится к обоснованию и разработке 
приемлемого для практического использования 
способа определения величины С ф . 


Математическая модель фонового 
химического состава поверхностных вод 


Рассмотрим одномерное уравнение турбулент- 
ной диффузии с функцией источника ДС)=-кС. 


дС_ + ѵ дС_ - д_ ( і) дС _ ч _ ^(2 
ді +ѵ дх ~дх [и дх ] кС ’ 


( 8 ) 


где х - координата пространства; ѵ - скорость дви- 
жения водных масс; И - коэффициент турбулент- 
ной диффузии. Пусть полное изменение расхода 
воды описывается ур. (9), а его частные изменения 
во времени и вдоль потока - выражениями (10, 11): 


сН т 

д О- 1() 


32 

дх 


( 9 ) 

( 10 ) 

(И) 


где г - время добегания водных масс; А, Я 2 — 
удельные скорости полного и частичного измене- 
ния расхода воды. Пусть А 2 /ѵ»соп8І. Тогда уравне- 
ние для полного дифференциала концентрации ве- 
щества можно записать в виде: 

сІС _ЗС , . ,^_дС^ ЗСдО._дС т \ 


дх 


д<2 д( 3(2 дх 3(2 

а выражение, описывающее диффузию, - в виде 
выражения: 

Г)д^ = І )( д_( 5СМ Л 


дх 2 


о(^) 2 2 


фх 3(2 дх 
г 


д 2 С п , дС 

ае 


(13) 


С учетом (12, 13) приведем ур. (8) к виду: 


лѴ , В(—) 2 - Я 

/о2 д С , /о оС ѵ 

да 2 0( я 2)2 


-^=о. 


(14) 


Полученное ур. ( 14) в зависимости от значений 
Д Я, Я 2 /ѵ, к имеет ряд известных решений. В част- 
ности, одно из решений при условии ДЯ»соп8І: и 
/иФ 0 принимает вид: 

С = С2 0 ' Н1 ~ Х+ ^ + С 2 д 0,н1 ~ х ~^, ( 15 ) 

где С, и С 2 - постоянные интегрирования; 


X - - 


Т>Дф 2 -А 

В(Щ 2 


В = (1-Я) 2 + 4- 


0(Щ 2 


Используя приближение \Іа 2 +х » а + уЦ 


ап- 


проксимируем ц выражением (16) и перепишем 
(15) в виде (17): 


В : 


_Я + 2к_ 

Ят \2 Я 


(16) 




с = сд 


+ С 2 ф 


(17) 


Принимая С[=0, можно получить частное реше- 
ние ур. (17) в виде: 


С = Со\ 


(18) 
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где С 0 и 0 О - начальные концентрация вещества и 
расход воды. 


Ур. (18) без специального обоснования исполь- 
зуется в гидрохимической практике при определе- 
нии фоновых концентраций веществ [6] и для дру- 
гих целей. На самом деле оно имеет вполне понят- 
ный физический смысл и описывает отклонение 
концентрации вещества от ее условно равновесно- 
го значения, соответствующего среднему геоме- 
трическому содержанию [7, 8]. Для подтверждения 
этого вывода рассмотрим на феноменологическом 
уровне систему «вода - порода», в которой проте- 
кают необратимые гидрохимические процессы. В 
целом, состояние этой системы неравновесно, но 
отдельные ее части квазиравновесны и имеют тер- 
модинамические параметры, время изменения ко- 
торых больше времени релаксации в этих частях и 
намного меньше времени, за которое устанавлива- 
ется равновесие во всей системе. В отдельной ма- 
лой подсистеме протекает I химических реакций, 
которые можно объединить в одну суммарную ре- 
акцию, описываемую выражением: 

Г ь \ 


А От 


КТ 


Іп П , - Іп Кг , 


(19) 


V 


где А Су и К? - общее изменение свободной энергии 
системы и суммарная константа равновесия при за- 
данной температуре Т; Д - суммарное произведение 
активностей компонентов, участвующих в каждой из 
реакций (Ь + 1 - количество веществ, участвующих в 
суммарной реакции). Используя упрощенное ура- 
внение Дебая-Хюккеля для описания связей между 
активностью иона а и его концентрацией С, после- 
довательно преобразуем (19) по схеме [8]: 

а=г\л[іС, 
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где г] - коэффициент, учитывающий заряд и эф- 
фективный диаметр иона; I — ионная сила раство- 
ра; Ъ - стехиометрический коэффициент; В , - кон- 
станта. При условии, что вероятность величины 
Іп С стремится к 1 /М, для математического ожида- 
ния логарифма концентрации рассматриваемого 
вещества Д1п С) будет справедливо: 


Е(\пС х ) = В,-^В ѵ Е(\пС ѵ ). (20) 

ѵ 

Таким образом, можно сделать вывод о том, что, 
во-первых, среднее геометрическое значение кон- 
центрации вещества соответствует условно равновес- 
ному состоянию системы «вода - порода», соответ- 
ствующему минимуму свободной энергии. Отклоне- 
ния от этого состояния определяются интенсивно- 
стью водообмена, регулирующей циркуляцию веще- 
ства и энергии в водном объекте в соответствии с за- 
висимостями вида (15, 17, 18). Во-вторых, под фоно- 
вым следует понимать условно равновесное состоя- 
ние водного объекта, сформировавшееся под влия- 
нием природных и антропогенных факторов за пред- 
шествующий статистически однородный период. 

Расчет нормативов ПДС 

С учетом уравнений (7, 20) последовательность 
расчета нормативов ПДС выглядит следующим об- 
разом: 

- по данным режимных гидрохимических иссле- 
дований (за весь период наблюдений) опреде- 
ляются средние геометрические значения и 
средние квадратические отклонения концен- 
траций нормируемых веществ; 

- с учетом различий гидрохимического режима 
различных типов водных объектов выявляются 
географические закономерности изменения 
средних геометрических значений и средних 
квадратических отклонений концентраций ве- 
ществ отдельно в водах больших, средних и ма- 
лых водотоков, озер, водохранилищ, прудов, бо- 
лот: для больших рек находятся зависимости из- 
менения содержаний веществ по длине водото- 
ка; для средних рек составляются карты (схемы) 
их изолиний; для прочих водных объектов — кар- 
ты (схемы) однородных гидрохимических райо- 
нов (обоснование различий гидрохимического 
режима больших, средних и малых рек на приме- 
ре бассейна Средней Оби приведено в [9]); 

- для проектируемых объектов при отсутствии 
необходимых наблюдений за составом сточных и 
поверхностных вод кратность разбавления прини- 
мается как п= 1 (величина ПДС корректируется в 
конце отчетного года по материалам гидрохимиче- 
ских наблюдений), для прочих объектов параметр 
п находится по формуле (4) по данным за предыду- 
щий период действия нормативов ПДС (до 5 лет); 

- для контрольных створов определяются фоно- 
вые содержания веществ и их средние квадрати- 
ческие отклонения: для больших рек - по зави- 
симостям от протяженности водотока или путем 
интерполяции между пунктами режимных на- 
блюдений; для средних рек - по картам (схемам) 
изолиний содержаний нормируемых веществ; 
для прочих водных объектов - по картам (схе- 
мам) однородных гидрохимических районов; 

- проводится расчет С дк ію ур. (7), а затем - вычи- 
сление величины ПДС. 
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Выводы 

Предложены новые способы расчета фоновых 
концентраций веществ в поверхностных водах и их 
нормативов ПДС в водные объекты со сточными 
водами. Преимущества данных способов по срав- 
нению с существующими заключаются в наличии 
теоретического обоснования, отсутствии жесткой 
привязки к ПДК, простоте расчетов и наглядности 
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представления материалов. Это позволяет более 
объективно определять допустимый уровень ан- 
тропогенного воздействия на реки, озера и болота. 
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